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CN . Abstract 

^ I On the basis of the decay couphngs /o — > vrvr, KK, rjrj, rjrj' , which had been found 

' before, in the study of analytical {IJ^^ = 00"'"^)-amplitude in the mass range 450- 



1900 MeV, we analyse the quark-gluonium content of resonances /o(980), /o(1300), 
C\i ! /o(1500), /o(1750) and the broad state /o(1420 t'j^^^)- The i^-matrix technique 

used in the analysis makes it possible to evaluate the quark-gluonium content both 
for the states with switched-off decay channels (bare states, /q"^*^) and the real 
60 . resonances. We observe significant change of the quark-gluonium composition in 

I the evolution from bare states to real resonances, that is due to the mixing of 

QQ ' states in the transitions fo{mi) real mesons /o("T'2) responsible for the decay 

. processes as well. For the /o(980), the analysis confirmed the dominance of qq 

component thus proving the nn/ss composition found in the study of the radiative 
decays. For the mesons /o(1300), /o(1500) and /o(1750), the hadronic decays do 
not allow one to determine uniquely the nn, ss and gluonium components providing 
relative pecentage only. The analysis shows that the broad state /o(1420 II^^q^) can 
mix with the flavour singlet qq component only, that is consistent with gluonium 
• origin of the broad resonance. 

> : 
•1—1 , 

^ : 1 Introduction 

The present paper continues the investigation of {IJ^^ = 00"'"+) resonances started in 
1^: we analyse the A'-matrix solution in which scalar glueball located near 1600 MeV 
appears, before mixing with neighbouring qq states due to decay prosesses. 

In [m, on the basis of experimental data of GAMS group [Q, Crystal Barrel Col- 
laboration 0] and BNL group 0, the if-matrix solution has been found for the waves 
00++, 10++, 02++, 12++ covering the mass range 450-1900 MeV. Also masses and total 
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widths of resonances have been determined for these waves. The following states have 
been seen in the scalar-isoscalar sector, 

00++ : /o(980), /o(1300), /o(1500), fo{U20t\'S), /o(1750) . (1) 

In 0, the resonances /o(1300) and /o(1750) are referred to as /o(1370) and /o(1710). 
The broad state /o(1420t 70°) is not included into the compilation the broad state is 
denoted in as /o(1530t 2^0) that represent the mean value for three solutions found in 
the K-matrix analysis; here we discuss the solution only where primary glueball is located 
near 1600 MeV; in this way, we use the mass of broad state found in this solution. 

For the scalar-isovector sector, the analysis ||1| points to the presence of the following 
resonances in the spectra: 

10++ : ao(980), ao(1520) . (2) 

In the compilation the state ao(1520) is denoted as ao(1450). 
As to tensor mesons, the following states are seen: 

12++ : a2(1320), 02(1660), 

02++ : /2(1270), /2(1525) . (3) 

Although in the analysis |I[ the CERN-Miinich group data ^ were not included into 
fiting procedure directy, these data had been fitted in previous papers p, ^, and it was 
specially controlled that solutions found in [Q] and |^, ^ for the vrvr channel are in a good 
agreement with each other. 

For the states shown in (1), (2) and (3), the i^-matrix poles and i^'-matrix couplings 
to channels vrvr, KK, rjrj, rji]', TrvrTTTr have been found in The i^'-matrix poles are not the 
amplitude poles, these latter corresponding to physical resonances, but when the decays 
are switched off, the resonance poles turn into the K-maXhx ones. In the states related to 
the i^-matrix poles there is no cloud of real mesons which are due to the decay processes, 
that was the reason to name them "bare states". The i^-matrix couplings for transitions 
bare state — >■ tttt, KK, rjr], rjrj' found in |]l| confirm the qq nonet classification of bare states 
suggested in (jTU], 



Still, the i^-matrix analysis [0 does not supply us with partial widths of the resonances 
directly. To determine couplings for the transitions resonance — > mesons, auxiliary 
calculations should be performed to find out residues of the amplitude poles. Calculations 
of the residues have been carried out in p| for scalar-isoscalar sector, that gave us the 
values of partial widths for the resonances /o(980), /o(1300), /o(1500), /o(1750) and broad 
state /o(1420 ^70^°) decaying into the channels 7171,71117171, KK, rjrj, rjrj'. 

In this paper we use the decay couplings for the reactions /o 7T7c, KK,T]r],T]T]' as 
a basis for the analysis of quark-gluonium content of scalar-isoscalar resonances /o(980). 
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/o(1300), /o(1500), /o(1750), /o(1420 l7(f°), supposing for these states a three-component 
structure: nn = [uu + dd)/\/2^ ss, gluonium. We demonstrate that hadronic decays do 
not determine the weight of all components providing the correlation between them only. 

The knowledge of coupling constants of bare states f^"'^^ tttt, KK,r]rj,rjrj' makes 
it possible to trace the evolution of the qq and gluonium components in /g-mesons by 
switching on/off the decay channels, thus establishing constraints for these components 
in a studied resonance. 

The paper is organized as follows. 

Section 2, being introductive, gives general picture of resonances in the scalar-isoscalar 
sector in the mass range up to 2300 MeV and presents qq classification of the OO"^"*" states. 
On the basis of this classification we analyse the K-matrix solution with a primary glueball 
located near 1600 MeV. Correspondingly, we explore parameters found for the broad state 
in this solution, /o(1420lVo°). 

In Section 3 we analyse quark-gluonium content of the resonances /o(980), /o(1300), 
/o(1500), /o(1750) and broad state /o(1420 1 70°) basing on the rules of quark combina- 



torics for the /o states 112 



In Section 4, using the iiT-matrix representation for the 00^^+ amplitude, we study 
the evolution of the quark-gluonium content of resonances /o(980), /o(1300), /o(1500), 
/o(1750), /o(1420 l\o°) by varying gradually the strength of the decay channels. 

Concluding, we state that our analysis based on the study of the hadronic decays only 
is not able to fix the gg/gluonium content of /o mesons unambiguously. For qualitative 
estimate of qq and gluonium components, one needs to incorporate additional information 
into analysis, such as partial widths of the /o-meson produced in 77 collisions, rates of 
radiative decays with the /o-meson production, or the /o production ratios in the decays 
of heavy mesons. 



2 Resonances in scalar-isoscalar sector 

The i^-matrix analysis |l| performed for the scalar-isoscalar sector at 450-1900 MeV 
allows us to reconstruct analytical form of the amplitude in the region shown in Fig. 1 by 
dashed line, since the threshold singularities of the OO"*"^ amplitude related to channels 
TTTT, TrTTTTTr, KK, 1^7] , 1^7]' are correctly taken into account. The amplitude poles which 
correspond to the resonances (1), broad state /o(1420 1 70°) included, are located just in 
the area where analytical structure of the amplitude 00"^"*" is restored. 

Below mass scale of the /^-matrix analysis |I| there is a pole related to the light 
(T-meson: in Fig. 1 its position, M = (430 — ^320) MeV, is shown in accordance with 
the results of the dispersion relation A^/D-analysis [|T^ (the mass region allowed by this 
analysis is also shown in Fig. 1). 
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The pole related to the light cr-meson, with the mass M ~ 450 MeV, has been obtained 
in a number of papers, see ^ for details. In the decay — > tt~^tt~^7t~ the cr-meson mass 
was found to be M = 480 ± 40 MeV [jl9[; corresponding pole location is also shown in 
Fig. 1 by crossed bars. 

Above the mass region of the K-matrix analysis, there are resonances /o(2030), /o(2100), 
/o(2340) |20|,|1. 



2.1 Classification of scalar bare states 



In [0, [TTjl, in terms of bare states, the quark-gluonium classification of scalar particles 
has been suggested. The results of the K-matrix analysis []I| support this classification. 

The bare state being a member of a qq nonet imposes rigid restrictions on the K- 
matrix parameters. The qq nonet of scalars consists of two scalar-isoscalar states, /g"'"^(l) 
and /q°^^(2), scalar-isovector meson Qq"'^'^ and scalar kaon Kq"''^^ . In the leading terms 
of the 1/A^-expansion p2[, the decays of these four states into two pseudoscalar mesons 
are determined by three parameters only, which are the common coupling constant g, 
suppression parameter A for strange quark production (in the limit of precise SU (3) flavour 
symmetry A = 1) and mixing angle Lp for nn = [uu + dd) / \/2 and ss components in /q"''^: 

nn cos if + ss sin ip . (4) 

The mixing angle defines the scalar-isoscalar nonet partners /o"'"^(l) and /q"''"^(2): 

y.(l)-^(2) = 90°. (5) 

The restrictions imposed on coupling constants allow one to fix unambigously the basic 
scalar nonet: 

l^Agg : /o""'(720 ± 100), a^'^"'=(960 ± 30), Js:o^'^""(1220+/i° ), /o^"''^(1260 ± 30) , (6) 



as well as mixing angle for /q ^(720) and /q '''^(1260): 



-70° 



- 5° 
10° 



(7) 



The nonet l^P^qq in the form of (P) has been suggested in DTof, where the i^'-matrix re- 
analysis of the Kt[ data |2^ has been carried out (bare states and their couplings for the 
OO'*'"^ and 10"*"^ waves have been found before, in |^). 

To establish the nonet of the first radial excitations, 2^Pf)qq, appeared to be more 
difficult problem. The K-matrix analysis [|l| gives two scalar-isoscalar states in the region 
1200-1650 MeV, their decay couplings satisfying the requirements imosed for glueball; 
these states are as follows: /o^"^^(1230tfo°) and /o^''^^(1600 ± 50). To resolve this dilemma 
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one needs the systematization of qq states on the (n, M^) plot (n being radial quantum 
number of the meson and M its mass): the systematization suggested in |^T| definitely 



proves /g"'"^(1600 ± 50) to be superfluous state on the qq trajectory. Correspondingly, 
/o^'^^^(1230l^3^o°) and /o*'^'-'^(1810 ± 30) should be the qq states. 

Below we present arguments based on the detailed consideration of the (n, M^) plot, 
while now let us discuss the variant which satisfies the constraints given by the qq trajec- 
tories. 

The nonet 2^Poqq looks as follows: 

2'Poqq: ft'i^'^^^-'io), /o''^'-^(1810 ± 30), a^^(1650 ± 50), K',''^%lS85t!ol) , .o^ 

(j)[f bare 1^1230)] = 40° ± 8° . ^ ' 

The ii'-matrix analysis [|l|, together with previous ones @, 0, enable to reveal in the scalar- 
isoscalar sector the bare state /g'^'"^(1600) which is an extra one for the nonet classification 
l^Poqq and 2^Poqq. At the same time the couplings of /q°^'^(1600) to the decay channels 
TTTT, KK, rjr], rjr]' obey the requirements imposed on the glueball decay. This gives the 
reason to consider this state as the lightest scalar glueball, 

0++ glueball : /o^"'''^(1600 ± 50) . (9) 

The lattice calculations are in an reasonable agreement with such value of the lightest 



glueball mass [|J 



After the onset of the decay channels the bare states are transformed into real reso- 
nances. For the scalar-isoscalar sector we observe the following transitions by switching 
on the decay channels: 

ft'{720) ± 100) ^ /o(980) , 
/o^«-(1230l\^°) /o(1300) , 

/o^'^'-'^(1260 ±30) /o(1500) , (10) 
/o^«-(1600 ±50) /o(1420tV^o°) , 
/o^«^^(1810±30) ^ /o(1750) . 

The evolution of bare states into real resonances is illustrated by Fig. 2: the shifts of am- 
plitude poles in the complex-M plane correspond to gradual onset of the decay channels. 
Technically it is done by replacing the phase spaces pa for a = tttt, tttttttt, KK, 7]7], rjr]' in 
the i^'-matrix amplitude as follows: Pa C,Pa, where the parameter ^ runs in the interval 
0<^<1. At^— >0 one has bare states while the limit ^ ^ 1 gives us the positions of 
real resonances. 
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2.2 Overlapping of /o-resonances in the mass region 1200—1700 
MeV: accumulation of widths of the qq states by the glueball 

The occurrence of the broad resonance is not at all an accidental phenomenon. It orig- 
inated due to a mixing of states in the decay processes, namely, transitions fo{mi) —>■ 
real mesons — *• fo{fn2)- These transitions result in a specific phenomenon, that is, when 
several resonances overlap one of them accumulates the widths of neighbouring resonances 
and transforms into the broad state. 

This phenomenon had been observed in [|, ^ for scalar-isoscalar states, and the follow- 
ing scheme has been sugg ested in H, the broad state /o(1420;Vo°) 

is the descendant 



of the pure glueball which being in the neighbourhood of qq states accumulated their 
widths and transformed into the mixture of gluonium and qq states. In this idea had 
been modelled for four resonances /o(1300), /o(1500), /o(1420tV^o°) and /o(1750), by using 
the language of the quark-antiquark and two-gluon states, qq and gg: the decay processes 
were considered to be transitions /o — > qq, gg, correspondingly, the same processes real- 
ized the mixing of the resonances. In this model, the gluonium component was dispersed 
mainly over three resonances, /o(1300), /o(1500), /o(1420t 70°), so every state is a mixture 
of qq and gg components, with roughly equal percentage of gluonium (about 30-40%). 

Accumulation of widths of overlapping resonances by one of them is a well-known effect 
in nuclear physics . In meson physics this phenomenon can play rather impor- 



tant role, in particular, for exotic states which are beyond the qq systematics. Indeed, 
being among qq resonances, the exotic state creates a group of overlapping resonances. 
The exotic state, which is not orthogonal to its neighbours, after accumulating the "ex- 
cess" of widths, turns into the broad one. This broad resonance should be accompanied 
by narrow states which are the descendants of states from which the widths have been 
taken off. In this way, the existence of a broad resonance accompanied by narrow ones 
may be a signature of the exotics. This possibility, in context of searching for exotic 



states, was discussed in 



The broad state may be one of the components which forms the confinement barrier: 
the broad states after accumulating the widths of neighbouring resonances play for these 
latter the role of locking states. Evaluation of the mean radii squared of the broad state 
/o(1420tV^o°) and its neighbouring resonances argues in favour of this idea, for the radius 
of /o(1420tVo°) is significantly larger than that for /o(980) and /o(1300) 0, |3| thus 
making it possible for /o(1420l 70°) to be locking state. 
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2.3 Systematics of the scalar- isoscalar qq states on the (n, M^) 
plot 

As is stressed above, the systematics of qq states on the (n, M^) plot argues that the 
broad state /o(142017q°) and its predecessor, /q"''^(1600 ± 50), are the states beyond qq 
classification. Following [^, we plot in Fig. 3a the (n, M^)-trajectories for /o, Gq and 
Kq states (the doubling of /o trajectories is due to two fiavour components, nn and ss). 
All trajectories are roughly linear, and they clearly represent the states with dominant qq 
component. It is seen that one of the states, either /o(1420l7o°) /o(1500), is superfiuous 
for qq systematics. Looking at the (n, M^)-trajectories for bare states (Fig. 3b), one can 
see that just /q"'''^(1600) does not fall on any linear qq trajectory. So it would be natural 
to conclude that /o*"'^"(1600) is an exotic state, i.e. the glueball. 

Relying on the arguments given by the systematics of qq states on the (n, M^) plot, 
we can state that classification represented by (8) and (9) as a basis for further analysis. 
Lattice calculations support the solution (8)-(9): calculations give values for the mass of 
the lightest glueball in the interval 1550-1750 MeV [0. 

3 Hadronic decays, rules of quark combinatorics for 
couplings and estimation of the quark-gluonium 
content of resonances 

In this Section, on the basis of the quark combinatorics for the decay coupling constants, 
we analyse quark-gluonium content of resonances /o(980), /o(1300), /o(1500), /o(1750) and 
/o(1420l7o°)- We would like to bring attention of the reader to ambiguities which are 
inherent in analyses studying hadronic decays of resonances only. 

3.1 Quark combinatorial relations for the decay couplings 

Within 1/A^ leading-order terms [Q, hadronic decays of meson resonances are determined 
by planar diagrams. An example of the process is shown in Fig. 4a: the quarks fiying 
away from the initial qq state produce in a soft way (i.e. at relatively large distances, 
r ~ Rconfinement) the ucw pair of light quarks {uu, dd, or ss) and turn into white hadrons, 
thus making it possible for initial quarks to leave the confinement domain. In the limit of 
fiavour SU(3) symmetry, the production of all quarks is equivalent; still, heavier weight 
of strange quark results in a suppression of the production probability of ss pair. So the 
following ratio of production probabilities takes place: 

uu : dd : ss = 1 : 1 : X , (11) 
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with A =0.4 — 0.8 for meson decays ||I2| , For similar decays, the couphng constants 



(see Fig. 4) differ by the coefficient which depends on the mixing angle of the nn and ss 
components of the initial meson and parameter A [Q. These coefficients are shown in 
Table 1 for the decays /o vrvr, KK, rjrj, rjr]', ri'rj' . 

The planar diagram for the transition gluehall two mesons is shown in Fig. 4b. 
Below, to illustrate the estimations, we consider gluonium as the two-gluon composite 



system: the large value of the soft-gluon mass, nigiuon ~ 700 — 1000 MeV supports 
the model, though, we should stress, the results are of more general meaning. 

The planar diagram of Fig. 4b represents a two-stage transition gluonium —^qq—^ 
two mesons, and the second stage is similar to the decay of qq meson. Because of that 
the relations between couplings in the transition gluonium two pseudoscalars are 
governed by the magnitudes given in Table 1 , with fixed values of mixing angle for the nn 
and ss components, qq = nn cos (fc + ss sin <^G; which have been formed in the process 
gluonium qq. The angle ipc is determined by the parameter Ag entering the first stage 
of the process gluonium — >• qq with relative probability uu : dd : ss = 1 : 1 : \g- Then 



cos^G = \/^^. (12) 

The relations between coupling constants for the decays gluonium tttt, KK, rjr], rjrj' , r(v( 
are given in Table 1 as well. In principle, Ag may differ from the suppression parameter 
inherent to transitions qq — > two mesons. However, it looks reasonable to use, as the first 
approximation, the same value of A for both stages of the gluonium decay, because the 
coefficients for the transitions qq two mesons change rather weakly with A belonging 
to the interval 0.4-0.8. 

The sums of couplings squared for the decay transitions of the qq meson and gluonium 
are of the same order, as it follows from the rules of the 1/A^-expansion ||2^ {N = Nc = 
Nf where and Nf are numbers of colours and flavours). Let us denote the sum of 
couphngs squared for transitions of Fig. 4a-type, qq — state J^^^sons, as g"^^ and the 
corresponding value for the gluonium decay, gluonium J^^^^ons, as gliuonium- The 
values g'^q and ggiuonium can be represented as discontinuities of the self-energy diagrams 
Figs. 4c, d, with cuttings shown by dashed lines: the cut block stand for the couplings 
shown in Figs. 4a,b. In terms of the 1/A^-expansion, the diagram Fig. 4c (and gg^) is of 
the order of 



7^- ~ G^- -Nf ~ (13) 

yqq ^ qq—meson—fqq J ' \ 



because Ggg-.^eson^w- ~ Likewise, for ggi^^^i^^ determined by diagrams of Fig. 

4d-type, one has: 

q^, . ~ . AT? ~ l-I . (14) 
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Table 1: 

Coupling constants given by quark combinatorics for /o-meson and glueball decaying into 
two pseudoscalar mesons in the leading terms of expansion; ip is the mixing angle for 

nn and ss states, /o = nn cos (p + sssinc/?, and is the mixing angle for rj — r]' mesons: 
T] = nn cos Q — ss sin and r]' = nn sin + ss cos 0, where cos ~ 0.8 and sin ~ 0.6. 



Channel 


The gg-meson decay 

couplmgs m the 
leading terms ol 1/N 
expansion 


The glueball decay 

1 • 

couplmgs 
in the leading terms 
of 1/A^ expansion 


Identity 
factor 

in phase 
space 


7r°7r° 


g cos</7 /v^ 


G /V2 + A 


1/2 


7r"'"7r~ 


g cos(p j\/2 


G /V2 + A 


1 


K+K- 


g{\/2snnp + \/~\cosip) / \/S 


gVV(2 + A) 


1 




5f(\/2sin(/7 + •\/Acos</?)/-\/8 


G ^A/(2 + A) 


1 


VV 


g {cos^ cos if/ -\/2+ 

-\/A sin if sin^ j 


G (cos2 + Asin2 0) /V2 + A 


1/2 


r]r]' 


g sin cos ^cos (f / \f2— 
\f\ sin(^j 


G COS0 sin0 (1 - A) /V2 + A 


1 


rj'rj' 


g{si\?Q cos V?/ a/2 + 

\f\ sin ip cos^ 0^ 


G (sin2 + Acos2 0) /V2 + A 


1/2 
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The coupling for the transition gluonium — > two — gluon state, Ggiuonium^gg, is of the 

order of Ggluonium-^gg ~ l/Nc. The estimates of COUphngS Gqg-rneson^qq and Ggluonium^gg 

are done by using basic self-energy diagrams for composite systems, which are of the order 
of unity: in cases under consideration such are the diagrams of Figs. 4e,f. 

The two-stage nature of the decay gluonium two pseudoscalars is a source of 
ambiguities in the determination of the quark-gluon content of mesons, if we restrain 
ourselves by hadronic decays only. The matter is that the qq meson with the quark content 
nn cos if + ss sin at 99 ~ ipc has the same correlations between coupling constants for 
the transitions /o nn, KK,r]r],r]r]' as those for the glueball. At 0.4 < Xg ^ 0.8 we 
have 24° < (fc ^ 32°: it means that the analysis of hadronic decays cannot distinguish 
between qq state with ip c:=npQ and true gluonium. 

3.2 Decay couplings for the resonances /o(980), /o(1300), /o(1500), 
/o(1750), /o(1420t^§0) into channels tttt , K K , r]r] , r]r]' , tttttttt 

The i^-matrix fit to the data gives us directly the characteristics of the bare states only. 
To extract resonance parameters one needs additional calculations to be carried out with 
the obtained amplitude. The couplings for resonance decay are extracted by calculating 
residues of the amplitude poles related to the resonances 0. In a more detail, the ampli- 
tude Aa^b, where a, b mark the channels nn, KK, rjr], rjrj', TrTrvrvr, can be written near the 
pole as: 

A., ^ 4^e^(^^"'+^^"n 5., . (15) 

The first term in ([T5| ) represents the pole singularity and the second one, Bab, is a smooth 
background. The pole position s = /i^ determines the mass of the resonance, with total 
width /i„ = Mn — iTn/2, and the real factors (7^"^ and g^^^ are the decay coupling constants 
of the resonance to channels a and b. The couplings g^^ given in Table 2 stand for the 
solution with glueball in the vicinity of 1600 MeV; also are shown the couplings for the 
predecessor bare states. 

The couplings of Table 2 demonstrate a strong change of couplings during the evolution 
from bare states to real resonances. Note that the change occurs not only in absolute 
values of couplings but in relative magnitudes as well. The resonances /o(1300), /o(1500), 
and /o(1750) demonstrate a reduction of relative weight of the coupling constant squared, 
9kk^ while the same coupling in the broad state /o(1420l7o°) increases. 

The growth of relative weight of gj^j^ in /o (1420^70°), that is the glueball descendant, 
can be unambigously interpreted: by accumulating the widths of neighbouring resonances 
this one acquires a noticeable qq component, with a large amount of ss state. 

Let us look at what the quark combinatorics rules given in Table 1 tell us about the 
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Table 2: 

Couplings squared, gl, (in GeV^ units) for bare states and their resonance-descendants 



State 


2 2 2 2 2 
9mr 9 KK 9rm 9nri' fl'TTTTTrTr 


2 


/o(980) 


0.167 0.528 0.06 - 
0.076 0.186 0.072 - 0.004 


0.755 
0.338 


/o^--(1220l^^) 
/o(1300) 


0.083 0.099 0.025 - 0.517 
0.026 0.002 0.003 - 0.132 


0.724 
0.163 


/o^"-(1265l^^) 
/o(1500) 


0.553 0.271 0.045 0.032 
0.038 0.009 0.007 0.006 0.074 


0.881 
0.134 


/o^'^''^(1820±40) 
/o(1750) 


0.059 0.019 0.001 0.043 0.262 
0.086 0.003 0.009 0.028 0.117 


0.384 
0.243 


/o(1420l^^°) 


0.124 0.062 0.025 0.008 0.924 
0.304 0.271 0.062 0.016 0.382 


0.513 
1.035 
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proportion of ss, nn and gluonium components in the studied resonances. The couphng 
constants squared for /o tttt, KK, rjrj, t]T]' can be written as follows: 




Sk/c = 2 I |(sm 'f + ^^ms<p)+G^ 

2 1 / , COS^ /-- . • 2 \ ^ I 2 i~\ \ • 2 \ \ 

9riv = 2 r ( ^ cosy? + VA sin(^ sm 9 ) + ^ (cos 6 + Asm-G ) I 

^w' = sin^ e cos^ e (^g {^cos(p - y/X sinip) + G ■^=j^ ■ 

The term proportional to g is responsible for the transition qq — > two mesons, while 
that proportional to G stands for the transition gluonium — > two mesons. Correspond- 
ingly, the magnitudes g"^ and G^ are proportional to probabilities to find the qq and 
gluonium components in the considered meson. 

First of all, let us determine the mean value of mixing angle for nn/ss components in 
the intermediate state, < (p >: 

fo ~^ gluonium + qq —>■ nn cos < ip > +sssin < ip >^ two mesons. (17) 

We define < p > a.s the angle for the coupling constants squared (|TBp at G = 0. Then we 
have for the studied resonances: 

/o(980) : <p>- -67° , (18) 

/o(1300) : < >~ -5° , 

/o(1500) : < ^ >~ 8°, 

/o(1420;^^°) : < ^ >~ 37° , 

/o(1750) : <p>- -27° . 

The value < v9[/o(1420l7o°)] >— 37° is very close to the mixing angle of the flavour 
singlet state, Psingiet = 35.3°. We see that /o(1420l7o°) is a mixture of the gluonium and 

singlet- 

Generally, by fitting formulae (|T^ to the coupling constants squared of Table 2, one 
can determine the mixing angle p as a function oiG/ g. The curves in Fig. 5 demonstrate 
the dynamics of p with respect to the ratio G/g for the resonances /o(980), /o(1300), 
/o(1500), and /o(1750). 
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The magnitudes g"^ and G"^ are proportional, correspondingly, to the probabilities for 
quark/gl uonium components, Wqq and Wgiuonium, to be in the considered resonance: 

9 ~ 9qq^qq i ^ ~ 9 gluonium^gluonium ■ (19) 

According to the rules of expansion p2|, the coupling constants, g^^ and ggi^ 



•Luoniumi 



are of the same order (see Section 3.1), therefore we accept as a rough estimation: 

G^/g'' = Wgluonium/Wqq . (20) 

In Fig. 5 we vary G/ g in the interval —0.8 < G/ g < 0.8; that, in accordance with (^), is 
related, to a possible admixture of the gluonium component up to 40%: Wgiuonium < 0.40. 

The qq components in the resonances /o(1300), /o(1500) reveal rather moderate change 
of versus a persentage of gluonium component: 

Wgi„m[h{l^m < 40% : -25° < <^[/o(1300)] < 13° , (21) 

and 

Wgiuomum[foil500)] < 40% : -3° < ¥^[/o(1500)] < 17° . (22) 
More sensitive to the glueball admixture is the qq components in /o(1750): 

Wgiuonium[fo{1750)] < 40% : -55° < </^[/o(1750)] < -2° . (23) 

The nn/ss ratio in /o(980) is also rather sensitive to the presence of gluonium component. 
However, for this case it is hardly possible to assume the glueball admixture to be more 
than 20% []TB|. Correspondingly, we have 

W^,,„„[/o(980)] < 20% : -95° < <^[/o(980)] < -40°. (24) 

The performed analysis clearly demonstrates the impossibility to fix unambigously the nn, 
ss and gluonium components in the resonances /o(980), /o(1300), /o(1500), /o(1750) by 
using their couplings to hadronic channels only. For more understanding of the structure 
of these mesons, one needs additional data. 

Still, the i^'-matrix analysis [|l| provides us with a bulk of information on meson states; 
that makes it possible to trace the evolution of coupling constants from bare state, /q"''^, 
to real resonances. 



4 Evolution of coupling constants for /o tttt, KK, 
TjiTj^ rjiTj' at the onset of the decay channels 

In this section, basing on the /^-matrix analysis results [|l|, we study the dynamics of pa- 
rameters for the resonances /o(980), /o(1300), /o(1500), /o(1750), /o(1420l^^°) by switch- 
ing on/off gradually the decay channels. 
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4.1 The onset of decay channels in the matrix ampUtude 

When the decay channels are switched off, the couphng constants to two-meson channels 
are determined by the residues of the i^-matrix poles, while, after switching them on, the 
coupling constants are determined by the residues of amplitude poles. Let us clarify this 
point in more detail. 

Scattering amplitude of the two pseudoscalar mesons had been fitted in in the form: 

A = k ^—^ , (25) 

I-ipK 

where K is the 5 x 5- matrix for five channels under investigation (vrvr, KK, 7777, 7777', mnni). 
The i^-matrix is real and symmetrical Kab{s) = Kf,a{s), I is the unit matrix, / = 
diag{l, 1, 1, 1), and p is the diagonal matrix for phase space, p = (im(7(p^^,p^/;-,/9^^,p^^/. 

The i^'-matrix element was represented in ||l[] as a sum of pole terms and smooth 
background fab{s): 

Kab{s)=Y.^\^ + fab{s) , (26) 
n f^n ^ 

Pn being the mass of bare state and g^^ the coupling of fo"'^^\pn) to the channel a. 

Fitting to data performed in fixes the parameters of the i^-matrix amplitude (cou- 
pling constants gab{s), masses of the A'-matrix poles pn and regular terms fat)- With 
the completely determined parameters we can investigate the dynamics of the onset of 
decaying processes. 



Such an investigation of the resonance evolution was suggested in |]35|, where rather 
simple variant had been considered: in the amplitude (|25|) the substitution p — > xp has 
been done, with x varying in the interval < x < 1. Thus constructed amplitude A{x) 
gives us the real amplitude at x = 1, and at x — the amplitude A(x) turns into the 
K-matrix one, A{x — *■ 0) — *• ii', which is the amplitude for bare states. Generally, the 
onset of the decay channels in the i^T-matrix amplitude can be investigated by substituting 
parameters as follows: 

U^)g^:^ , fab ^f{x)fab . (27) 

Here the parameter-functions ^n{x), C/(x) obey the requirements ^„(0) = ^/(O) = and 
(^„(1) = ^/(l) = 1. A simple variant studied in ||35| corresponds to C,n{x) = a/x and 
^/(x) = X. 

The parameters C,n{x) and ^/(x) control the dynamics of the onset of the decay channels 
for resonances. This dynamics may be different for different states, say, qq state and 
gluonium. Below we use ^„(x) = a/x for the states connected with /o(980), /o(1300). 
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/o(1500), /o(1750), while for the broad state the dependence Cn{x) = x^l'^ is accepted. 
For the background term we use ^n{x^ = x. 

4.2 Evolution of the decay couplings 

Below the evolution of states related to the resonances /o(980), /o(1300), /o(1500), /o(1750), 
and /o(1420l^^°) will be considered one by one. The procedure is as follows: the substi- 
tution (^) is done and x takes the numbers x = 0.1, 0.2, 0.8, 0.9. The value x = 
corresponds to bare states, and the couplings and masses had been found as fitting param- 
eters in 1^ ; for x = 1 calculations were performed in and the couplings are presented in 
Table 2. Furthermore, for different but fixed x we find the position of pole and calculate 
the residues of the amplitudes thus determining g^^\x). 

4.2.1 Resonance /o(980) 

For the /o(980) the normalized couplings 7a = Qa/ \J qI-,, + 9\:ki where a = 7nr,KK are 
shown in Fig. 6a. Correspondent poles at different x are placed on the trajectory 
/q°^'^(720) — /o(980) (see Fig. 2). Looking at Fig. 6a, one can see that at small x 
gl^ < g]^j^, that is natural, for the state /q"'"^(720) is close to the flavour octet. With 
the increase of a;, the coupling constants become equal to each other, and in the in- 
terval 0.6 < X < 0.8 the coupling to pion channel is greater than to kaon channel, 
qIt^^x ~ 0.7) > gl^j^ix ~ 0.7), thus revealing a relative reduction of ss component. How- 
ever, at X ~ 0.8, when the amplitude pole approaches the KK threshold, see Fig. 2, the 
relative weight of the KK channel is strenthened to certain extent. 

For every fixed x, the formula (16) has been fitted to the values of gl^^ix) and gl^^^^) 
in order to find out as a function of G/g: a set of obtained curves is shown in Fig. 6b. 
At X = the curve ip{G/g) corresponds to maximal possible values of mixing angles; with 
the increase of x these curves move smoothly to smaller \(p\, while at x > 0.8 they turn 
backward rather sharply. The curve for to x = 1 is also shown in Fig. 5, it coincides with 
a good accuracy with that at x = 0. 

In Fig. 6b there are also two curves for if with gradual accumulation of the gluonium 
component which results in Wgiuonium = 20% at x = 1: recall, we suggest Wgiuonium < 20%. 
These are dotted curves which originate from the point G/g = at x = (in the K-matrix 
analysis [|] the bare state fl^"'^ (720) is the pure qq) but furthermore they drift to G/g > 
and G/g < and end at x = 1 with G/g = 0.45 and G/g = —0.45, correspondingly. 
These final points stand for Wgiuonium = 20%. Respectively, the variations of ip are shown 
separately by Fig. 6c; the area inbetween the curves is just the region of reasonable values 
of (fi with the evolution of the decay widths. The lower curve tells us that the reduction 
of gl^j^ with the increase of x to the region x < 0.8 is plausible due to the growth of the 
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gluonium component leaving the ratio nn/ss approximately constant. The upper curve 
related to G/g < testifies the possible evolution when, in parallel with the increase 
of the gluonium component, the weight of the ss component becomes smaller. But at 
X > 0.8, both curves demonstrate sharp increase of the ss component. 

At X = 1, when mixing angle ip which defines the quark content of the /o(980), 
(nncos + ss sin 0), its value may vary in the intervals —90° < (p < —48° and 85° < f < 
90° (Fig. 6c). It would be instructive to compare these values with mixing angles obtained 
for radiative decays 0(1020) — > 7/o(980) and /o(980) — > 77. The combined analysis [0 
for these decays provided us with two possible solutions: cp = —48° ±6° and = 86° ±3°. 

The areas for y)[/o(980)] allowed by radiative and hadronic decays are shown in Fig. 
7. One can see that the constraints for mixing angle obtained from the study of hadronic 
decays of the /o(980) are in a nice agreement with the values obtained from the study 
of radiative decays. However, one should stress, that we cannot expect from hadronic 
processes more rigid limitations for mixing angle of the /o(980) than those known from 
radiative decays. 



4.2.2 Resonances /o(1300), /o(1500) and /o(1750) 

For the resonances /o(1300), /o(1500) and /o(1750), the fraction of the ss component 
decreases. They fiow away from these resonances and enter the broad state /o(1420l7o°). 

Figure 8 demonstrates normalized coupling constants squared 7a = Qa/ \Jj2b9b ^ func- 
tion of X. One can see that for /o(1300), /o(1500) the values 7^^ are falling down while 
77r7r increase. At the same time, for the broad state the normalized coupling jj^j^ is 
growing up. 

The description of ratios 7a by Eq. (16) is shown in Fig. 8: one can see that quark 
combinatorial rules describe reasonably a number of data on the decay coupling constants. 



But, as was stressed above, Eq. (|T6D does not define the content of a resonance, providing a 
correlation only between mixing angle p and ratio G/g. These correlations are presented in 
Figs. 9a,b,c for the /o(1300), /o(1500), /o(1750) at various x. Implying the predecessors of 
these resonances to be pure qq states, we display the correlations {p, G/g) at G'^/g'^ < 0.4 



that corresponds to (|21|), (p^ and (p^). Dotted curves in Figs. 9a,b,c stand for {p,G/g) 



related to a maximal capture of the gluonium component by these resonances. 



4.2.3 Broad state /o(1420l^^°) 

The evolution from the /o^'"^'^''(1600) to /o(1420I'^7q°) is accompanied by the accumulation 
of qq component and the growth of ratio g/G. In this way. Fig. 9c demonstrates the 
correlation {p, G/g): for small x, when parameters of the corresponding state are close to 
those of /o^"'^'^''(1600), the ratio g/G is small. We see that for the broad state the mixing 
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angle (y? does not depend practically on gjG at jC^ < 0.50. 

The value of the mixing angle at x = 1, Lp = 37°, proves that the /o(1420t7o°) 
can accumulate the flavour singlet component of qq only; that perfectly agrees with its 
gluonium origin. 

Figure 10 demonstrates the change of (f with increasing x: this curve does not depend 
on the rate of accumulation of the (qq) singlet component. 

5 Conclusion 

We have performed the analysis of coupling constants for the resonances /o(980), /o(1300), 
/o(1500), /o(1750), and the broad state /o(1420t7o°) to channels mr, KK,rir],r]r]' as well 
as observed the evolution of bare states into these resonances by switching on/off the 
decay channels (see (10) and Fig. 2). Our analysis has been based on Solution 11-2 of 
the paper [0; in this solution, the bare state /o^"^'^'*(1600 ± 50) is the candidate for the 
glueball. 

During the evolution of states, the coupling constants /o — * hadrons change consid- 
erably not only in magnitude but relative weight as well, that is due to a strong mixing 
of states, because of the decay processes /o(w^i) real mesons /o(^2)- Using the 
language of bare states , fo°'^'^\ a rigid classification can be established for qq states; 
however, for real resonances the presence of the gluonium component results in certain 
uncertainties. 

For the discussed resonances the results are as follows. 

1. /o(980): This resonance is dominantly the qq state, the admixture of the glueball 
component is not more than 20%, Wgiuonium < 0.20. Rather large ss component is also 
present. Taking account of the representation qq = nncosip + sssiny^, the hadronic 
decays give us the following constraints: —90° < < —40° or 85° < (p < 90°, which are 
in agreement with data on radiative decays /o(980) 77 and 0(1020) — >• 7/o(980) 0] 
(see Fig. 7). Rather large uncertainties in the determination of mixing angle are due to the 
sensitivity of coupling constants to plausible small admixtures of the gluonium component. 
When the gluonium component is absent, hadronic decays provide ip =< ip >= —67°. 

2. /o(1300) and /o(1500): These resonances are the descendants of bare qq states 
^^(bare) ^^230+150) ^nd /q^^'"^^^ (1260 ± 30) which both of them are flavour singlets. The 
resonances /o(1300) and /o(1500) are formed due to a strong mixing with gluonium state 
/q"'^'^(1600) as well as with one another. The qq content, qq = nncosip + sssin(y9, in 
both resonances strongly depends on the admixture of the gluonium component. At 
Wgiuonium < 0.40 the mixing angles change, depending on Wgiuonium, in the intervals 
-25° < (/?[/o(1300)] < 13° and -3° < (/?[/o(1500)] < 17°. 

3. /o(1750): This resonance is the descendant of bare state 2^Piqq, /q^''"''^^ (1810 ±30), 
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and this bare state has the flavour wave function close to the octet one. During evolu- 
tion and mixing (presumably with the gluonium) the quark component, qq = nncosy? + 
sssiny), can change significantly: with Wgiuonium < 0.2, we have —45° < (p[fo{1750)] < 
—10°. Therefore, /o(1750) keeps his large ss component, and rather small coupling con- 
stant /o(1750) — i> KK should not mislead us, for the production of KK is suppressed at 

ifi 30°, see Table 1 and Eq. (|l|). 

4. /o(1420l^g°): This broad state is the descendant of the /o^''"'''^(1600 ± 50) which 
we believe to be the glueball. The analysis of hadronic decays of this resonance confirms 
the glueball nature of this resonance: the qq component is allowed to be in the flavour 
singlet only, (^[/o(1420I'^7q°)] ~ 37°, though the value of the possible admixture of {qq)singiet 
cannot be fixed by hadronic decays. In this way, the /o(1420l7o°) is a mixture gluonium + 
(qq) singlet] the impossibility to find out the quark-antiquark component is due to the fact 
that correlations between the decay coupling constants are the same for the gluonium and 
qq singlet. 
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Figure 1: 

Complex M-plane for the {IJ^*^ — 00+"*") mesons. Dashed hne encircles the part of 
the plane where the X-matrix analysis [1] reconstructs analytical X-matrix amplitude: 
in this area the poles corresponding to resonances /o(980), /o(1300), /o(1500), /o(1750) 
and the broad state /o(1420 ^jq^) are located. Beyond this area, in the low-mass re- 
gion, the pole of the light cr-meson is located (shown by the point the position of pole, 
M — (430 — i320) MeV, corresponds to the result of N/D analysis [11]; the crossed 
bars stand for cr-meson pole found in [19]). In the high-mass region, one has resonances 
/o(2030),/o(2100),/o(2340) [17,18]. Solid lines statnd for the cuts related to the thresh- 
olds TTTT, TTTTTTTr, KK, Tjr], 7)7)'. 
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Figure 2: 

Complex M-plane: trajectories of the poles for /o(980), /o(1300), /o(1500), /o(1750), 
/o(1420 lyo^^) during gradual onset of the decay processes. 
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Figure 3: 

Linear trajectories in (n, M^)-plane for scalar resonances (a) and scalar bare states 
Open points stand for predicted states. 
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e f 
Figure 4: 

(a,b) Diagrams responsible for the transition of scalar-isoscalar state into two pseudoscalar 
mesons, /o PP, in the leading 1 /A^^c-expansion terms. (c,d) Self-energy diagrams which 
determine dl i^oi quarkonium) and J2 (for gluonium); the cutting is shown by dashed 
line. (e,f ) Self-energy diagrams which determine the order of value of couplings for the 
transitions quarkonium — > qq and gluonium — > gg. 
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Figure 5: 

{G/g,(p)-plot: the trajectories show correlations between G/g and (p for /o(980) (sohd 
hne), /o(1300) (long-dashed line), /o(1500) (dased line), /o(1750) (dotted line). 
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Figure 6: 

The evolution of the /o(980) resonance parameters by switching on/off the decay channels 
(x = corresponds to the bare state while x = 1 stands for the real meson), (a) Nor- 
malized coupling constants g-,,.,,/ \J gl^, + 9kk 9kk/ \J qI-k ~^ 9kk ■ (^) '^^^ correlation 
curves {G/g, if) at ten fixed values of ,t; dotted curves stand for maximal accumulation of 
the gluonium component, Wgiuonium = 20%. (c) The evolution of ip in the qq component 
qq — nn cos </? + ss sin (/? at maximal accumulation of the gluonium component: solid curve 
stands for positive G/g and dotted one to negative G/g. 
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Figure 7: 

The (</?, ii!jg(98Q))-plot: the shaded areas are the allowed ones for the reactions 0(1020) 
7/o(980), /o(980) ^ 77 and /o(980) ^ tttt, KK. 



27 



7a 

0.8 










fo(1300) 






7a 

0.8 








fo(1500) 


TT7T 




0.6 

0.4 


o- 

■ 


x> 
m 


"m - 


iS 

\ 

■ 


\ 


TTTT 

m 


- 

-■ 


0.6 
0.4 


• — •- 

0- ©■ 


* 

-e 




m m 9 

0- -0 -0 


KK 
a 

rij} 


- 

-s 


0.2 










-e 


-o-o- 





0.2 










□ 
















KK 


















































1 







0.2 




0.4 0.6 


0.8 


1 

X 


1 





0.2 


0.4 


0.6 


0.8 


1 

X 


0.8 










fo(1420_T 


) 

nn 




0.8 




• 




— • — • - 


fo(1750) 
• — • * 


7T7T 




0.6 


o 


o 


o 




O 0_ O _D 


KK 


p 


0.6 


"° n 








rjr]' 




0.4 
0.2 


• - 

-□ 


■- 

□ 


-■- - 

n 


■ 




VV 
~ ■" ■ 

7777' 


■ 


0.4 
0.2 


<R 

\ 

s 


-0 


-0- G- 


vv 


□ 
■ 

















□ 













KK 











0.2 




0.4 0.6 


0.8 


1 

X 







0.2 


0.4 


0.6 


0.8 


1 

X 



Figure 8: 

Evolution of the normalized couplings 'jada/ \l^b 9b within onset of the decay channels 
for /o(1300), /o(1500), /o(1420 tii^), /o(1750). Curves demonstrate the description of 
couplings by Eqs. (16). 
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Figure 9: 

Curves of correlations {G/g,^) for /o(1300), /o(1500), /o(1750) and {g/G,ip) for 
/o(1420 ^70^°). at ten fixed values of x; dotted curves stand for maximal accumulation of 
the gluonium component, Wgiuonium — 40%. 
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Figure 10: 

The evolution of </? in the qq component qq — cos </? + ss sin (p at maximal accumulation 
of the gluonium component, Wgiuonium = 40%; solid curves stand for positive G/g and 
dotted ones to negative G/g. For /o(1420 lyo^^), the if does not depend on the value of 
the adopted qq component. 
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